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Abstract. - We consider the inverse Landau-Zener problem which consists in finding the 
energy-sweep functions W(t) = £i(£) — £2(t) resulting in the required time dependences of 
the level populations for a two-level system crossing the resonance one or more times during the 
sweep. We find sweep functions of particular forms that let manipulate the system in a required 
way, including complete switching from the state 1 to the state 2 and preparing the system at 
the exact ground and excited states at resonance. 

The Landau-Zener (LZ) effect jl[ ||, ||, which consists in the quantum-mechanical transition 
between the two time-dependent levels of a two-level system (states 1 and 2) caused by crossing 
the resonance, is a well known phenomenon in many areas of physics. Well known applications 
of the LZ effect are those to molecular collision and dissociation^, Recently, LZ effect has 
been used to measure the tunnel splitting in molecular magnets whereby detecting topological 
interference effects ^ which were predicted theoretically [Q, |], |j] ■ If the system was in state 

1 before crossing the resonance, then, for a constant sweep rate v, the probability to stay in 
this state after crossing the resonance is given by 

where v = \W{t)\, W(t) = e 1 (t)-e 2 {t), e lt2 (t) = (V>i, 2 |#lV'i,2>, A = 2{ifo.\H\fo), W{t) satisfies 
H^(±oo) = ±oo and A is the tunnel level splitting. For fast energy sweep rates v the system 
spends too little time in the vicinity of the resonance so that the tunneling matrix element 
(ipi\H\ip2) cannot bring the system into the state 2 (say, onto the other side of a potential 
barrier), and the system remains in the initial state 1. For slow sweep rates the probability to 
remain in the initial state becomes exponentially small, and the system travels into the state 

2 remaining on the lower exact energy term e_ (t) of 

1 r 



e±(t) = - [e 2 (t) + ei(t) ± y/W*(t) + A* 



(2) 

The corresponding adiabatic time dependence of the probability to be in state 1 is given by 



NWI'-l(*.l*-W)l» = ;|i- ^' +A , )- W 
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In the case where fiv ~ A 2 and v is constant, the solution of the Schrodinger equation for 
|ci(i)| 2 can be found analytically Jl], ||, [| and it is oscillating at t > slowly approaching P 
given by Eq. ([!]) at t — > oo. 

The known results above are valid for the linear energy sweep, W(t) = vt. One can consider, 
however, other sweep functions W(t) that behave nonlinearly in the vicinity of the resonance, 
W ~ A. If the tunneling resonance is not too narrow and its position is well defined, nonlinear 
sweeps of a controlled form can be realized experimentally. Although no general analytical 
solution of the Schrodinger equation is known for the LZ problem with an arbitrary W(t), one 
can look for the functions W(t) that result in a desired time dependence of the probabillity 
|ci(i)| 2 , thus solving the inverse LZ problem. It is the main motivation of this Letter to study 
this inverse problem which, for instance, will allow to determine W(t) needed to manipulate 
and to prepare a quantum system in a specific state. 

One class of solutions of the time-dependent Schrodinger equation are smooth, non- 
oscillating functions of time that describe the complete switching of the system from the state 
1 to the state 2, e.g., Eq. (||). Another class of solutions are those starting in ipi at t = — oo off 
resonance and ending in a given superposition of the exact states tp± at resonance [W(t) = 
for t > 0]. The most interesting variants of the latter are the system in the ground state ip- 
at resonance and in the excited state ip+ at resonance. These possibilities to manipulate the 
system in a required way could prove to be useful for the quantum computing (QC). Another 
interesting possibility is to sweep at a high rate (nonlinearly in time) far from the resonance 
thus making the whole process of crossing the resonance much shorter. Exploring these new 
types of solution of the LZ problem is the aim of this Letter. 

To this end, we consider the Schrodinger equation for the coefficients in the wave function 
ip{t) = ci(t)ipi + C2{t)ip2 that can be written in the form 

hp{t) = -~Aq{t), p(-oo) = l 

hq{t) = iW{t)q{t) - ~Ap(t), q(-oo) = 0, W(-Oo) = -oo (4) 

with p(t) = e 1£lt ci(t) and q(t) = e leit C2(t). We are looking for sweep functions W(t) that 
result in a required time dependence of the state populations. Before solving this inverse LZ 
problem analytically, we will illustrate numerically the possibility of time-symmetric solutions 
of Eq. (Q) going from p{— oo) = 1 to p(oo) = for strongly nonlinear sweep functions W(t). 
The simplest of these cases is the cubic-parabola sweep with W(t) = —2.5t + 1.83< 3 where we 
have set A = H = 1 in Eq. (0). The time dependence of the probability |p(i)| 2 shown in Fig. 
can be fitted by |p| 2 (i) = [1 — tanh(sinh(£))]/2 which is, however, not the analytical solution 
ofEq. (|). 

To find analytical solutions of the inverse LZ problem, one should rewrite Eq. (^) in terms 
of the probability m = \p\ 2 and some phase variable. We thus define 

p = \/mexp(iipp), q = vl — mexp(j<p 9 ) (5) 

and substitute it in Eq. (Q). Evidently the global phase of the system (ip p + f q )/2 is irrelevant 
and the resulting Schrodinger equation can be formulated in terms of m and the phase 
difference <p — ip q — ip p : 



Tim — A y/ m(l — m) sin <p 
A 1 - 2m 
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Fig. 1. - The cubic-parabola energy sweep with Wit) = —2.5t + 1.83t 3 and A = U = 1. Time 
dependence of the probability |ci(t)| 2 = \p(t)\ 2 can be fitted with |p| 2 = [f — tanh(sinh t)]/2 with 
graphic accuracy. The function \p\ = [1 — tanht]/2 is shown for comparison. 



This system of equations can be brought into a more elegant form by introducing 

x = 2m-l, -1 < x < 1, (7) 



y = arcsinx, —tt/2 < y < n/2, (8) 

and the reduced variables 

T = (A/h)t, W = W/A. (9) 
Then Eqs. (^J) can be put into the reduced form 

d T y — ship 

d T <p — W — tany cos tp. (10) 

The boundary condition for this equation is y{— oo) = 7r/2 whereas the value of ip(— oo) 
is irrelevant. In general, far from the resonance ip(r) strongly oscillates with time and the 
derivative d T y is not small. Below we will consider special kinds of solutions of Eq. (jl(^) which 
are characterized by a smooth, non-oscillating dependence sinter) and d T y{—oo) = 0. We will 
solve the inverse problem for Eq. ( |To| ) which consists in finding W(t) resulting in a required 
function x(t). One easily obtains from Eq. (|l(|) a general formula 

W(t) = d T ip + tan ycos tp = d^y /cos ip + tany cos 95, cos (p = ±\/l — {d T y) 2 . (11) 

The sign in front of cost/? determines the sign of the function W(t). Since |ci(— oo)| 2 = 1 and 
thus x{— 00) = 1 and y{— 00) = 7r/2 one has tany(— 00) = 00. To comply with the condition 
W{— 00) = —00 one has to choose, in general, the negative sign: 

W(t) = ~8 2 T y /y/l - (d T y) 2 - tany^l - (dry) 2 - (12) 

Before considering particular cases, we will comment on some general properties of the 
functions x(t) and W(t). If x(t) is an odd function of time then y(r) is also odd and d T y(r) 
is even, whereas d 2 y(r) and tany(r) are odd. Then from the first of Eqs. ( |l0| ) follows that 
sin^(r) is even. If the sign of cos(p does not change, which is the generic situation, cosp(r) 
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is even. Then from Eq. ( |IT| ) follows that W{t) is odd. In a similar way one can show that 
if x(t) is even (the system returns into the initial state ipi after crossing the resonance) then 
W(t) is even, too. 

There is an important special case, however, which is realized if the derivative d T y at some 
point to attains its maximal value \d T y\ — 1 which follows from the first of Eqs. (|To|) . At 
this point cos if turns to zero and changes its sign to positive. The latter corresponds to the 
phase ip going from the interval ir/2 < <p < ir into the interval < ip < tt/2. Note that 
d 2 y — at t — tq thus Eq. (|l^) does not diverge. If r = then cos <p(r) is an odd function: 
cos (p < for r < and cosp > for r > 0. As a result, the parity of W(t) becomes inverse 
to the parity of x(t). If t =/= then for a symmetric function x(t) the similar happens at 
r = —tq. That is, cos tp{r) changes sign two times and is an even function, thus the parity of 
W(t) coincides with that of x(t). It should be noted that parity considerations for the direct 
LZ problem do not apply: Symmetric forms of W(t) do not necessarily result in symmetric 
solutions x(t). If, however, the solution x(t) satisfies x(±oo) = =pl or ir(±oo) = ±1 and W(t) 
is time symmetric, then x(t) is also symmetric. This follows from the time reversibility of the 
Schrodinger equation and the symmetry of the boundary values of x(t). 

Another implication of this analysis is the existence of the maximal rate of the probability 
change. We have seen that \d T y\ < 1, then from d T x = \/l — x 2 d T y follows \d T x\ < 1. This 
limitation is physically clear: Since tunneling between ip\ and "02 is caused by the matrix 
element A = 2{ij)\\H \4>2), the latter sets up the natural scale for the tunneling rate, see Eq. 
(^|). To completely tunnel from ipi to "02, the system should spend some minimal time in the 
vicinity of the resonance, that is, the sweep cannot be too fast. Below we will discuss two most 
important families of the functions x(t) and corresponding W(t) which depend on r in the 
combination u — ar, the parameter a being the measure of the sweep rate. We will see that 
there is a maximal value a = et max for which \d T y\ — 1 at some r = tq. We will also see that 
although the functions x(t) do formally exist for a > a max , one cannot extend W(r) into this 
region, as a result of the physical limitation on the tunneling rate. Thus solutions x(t) with 
a > amax are not realizable. 

One of possible solutions is of the form 

u n 

x =(_!)" ii = dT . (13) 

For n odd it describes switching of the probability from ipi to ip2, whereas for n even the 
system returns into the initial state. For n = 1 the form of Eq. (^3|) corresponds to Eq. (||) by 
taking into account Eq. (0) with m — \c\\ 2 . However, the reader should note that Eq. j§) is 
realized with the very slow sweep W{r) of an arbitrary functional form, whereas Eq. (|l3) is 
valid for a particular form of W(j) which is not necessarily slow. One obtains from Eqs. (|^) 
and d||) 

1 runny 1- * 

d ^y = u o d r x = (- 1 ) 



2/j 



VT^x 1 i + u 

„ a 2 nu n ~ 2 \n-l-(n+ l)u 2n ] 
dly = (-1 n [ - 14 



whereas tany = x/yl — x 2 = (—l) n u n . The derivative d T y has a chance to attain the value 1 
at its maximum that is determined by d^.y = 0, i.e., at 

1 \ 1/(2") 

n — 1 x 



= ±1 — ) ■ (15) 
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At this point one has 



\{drV)r 




thus the condition 



a < a r 



2 

n+1 



(n-l)/(2n) 



2n/(«-l) 



n = 1 
n > 1, 

n = 1 
n > 1 



should be fulfilled. For a < a max 
(-1)"-^" 



W(r) = 



-MY 

1 \n+l J 

Eq. (p|) yields 

a 2 n [(2n + l)u 2 "- 2 - (n - l)u~ 2 ] 



(16) 



(17) 



yd - [ami"- 1 / (1 + u 2 ™)] 2 



(l + u 2 ™)' 



(18) 

This function is plotted for n = 1 and n = 5 and different values of the parameter a in 
Figs. |2|a, and For a — > a max discontinuities in Eq. (|l8|) are formed at it = ±it max . These 
discontinuities result from the invariably negative choice for cos (p. For a = a max the sign of 
cos (p can change at u = w max and then these discontinuities disappear. This form of W(t) 



includes the factor sgn (l + u 2n ) — (anu 11 x ) , in addition to Eq. (|18|). Note that Eq. (|l 

cannot be extended into the region a > a max preserving W(t) real, in contrast to Eq. (|l3| 
The small- and large-argument forms of W(t) are 



W{t) oc (-1)"" 1 
In the case n—1, Eq. (fL8|) simplifies to 



,71-2 



\u\ > 1 

H < i, 



n > 1 



1 - 3a 2 -I- 2u 2 



1 + tt 2 



2u 2 



(19) 



(20) 



One can see that for «»1 the function W(t) describes a linear energy sweep with the velocity 
v = a. For u ~ 1 the form of VF(r) is nonlinear and for a < a max = 1 the resonance is crossed 
three times, including u = 0. For a = 1 Eq. (Eoh yields the even function 



W(t) = 



2-2r 2 



(l + r 2 )V2+^ 



(21) 



This sweep crosses the resonance two times, in the forward and backward directions, which 
means that the system is in the excited state ip2 at the end of the sweep at t = oo. For a < 1 
and far from the resonance Eq. (^0|) has the form 

(22) 



W(t) = AW(t) = Aar = (A 2 a/Ti) t, 



thus the corresponding linear sweep rate is v = A 2 a/h. The maximal sweep rate that ensures 
the complete switching to ip2 is realized at a = 1 is u max = A 2 /h. For a comparison, in the 
regular LZ effect with the probability to stay in the initial state given by Eq. ([!]) the sweep 
rate v <C irA 2 / (2h) is required for a complete switching. For v — A 2 /% the probability to stay 
in the initial state is still P = 0.208. 

Another solution that describes switching from state 1 to state 2 is 



tanhu, 



<JT. 



(23) 
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Fig. 2. - a - Energy sweep function W(r) of Eq. (|lf]) with n = 1 [or of Eq. (53)] for different values 
of the sweep rate a. b - The same with n = 5. c - The population of the exact excited level, |c+| 2 , 
for n = 1 and different values of a. d - The same for n = 5. Dashed lines in all graphs correspond to 
a = a max and W(t) without jumps [see, e.g., Eq. (gj)]. 



Here the sweep function determined from Eq. (|12|) reads 

sinhu 



W(r) 



\J\ — a 2 J cosh 2 u 



1 



2a 2 



cosh 2 it 



(24) 



In the range l/y/2 < a < 1 the function W(t) changes its sign three times and its derivative 
at r = is negative. With a — -> 1 this derivative increases unlimitedly and a jump of W(r) at 
r = is formed. For a = a max = 1 the solution with eliminated jump is an even function 



W(t) = 



coshi 



cosh r. 



(25) 



This sweep function crosses the resonance two times. Again, Eq. (24) cannot be continued 
into the region a > 1. 

Now we return to the solution for the probability x described by Eq. (IF}!) and analyze its 
behavior for large n. It has a form of three plateaus x = 1 for u < — 1, x = for — 1 < u < 1, 
and x = — 1 for 1 < u with two narrow transition regions between them. In the central region 
— 1 < u< 1 one has | C2 1 2 — |ci| 2 = (1 + x)/2 = 1/2, i.e., the system is in the mixture of the 
states ipi and ip2 with equal weights. The sweep function W(t) of Eq. ( |la ) is fast approaching 
the resonance at u = ~~ 1 and remaining in the close vicinity of the latter until u = 1. Since 
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the probabilities |ci| 2 = | C2 1 2 = 1/2 do not oscillate in the range —1<u<1 (unlike the case 
in which the system is instantaneously brought into resonance), one is lead to the conclusion 
that the system is in one of the exact quantum states ip- or ip+. It is naturally of interest 
to question in which of these exact states does the system come under the influence of the 
energy sweep described by Eq. (|l^). We will see below that both variants can be realized. 
This opens a possibility to prepare the system in the ground or excited state in resonance by 
the appropriate choice of W(t). Indeed, for n > 3 one has W(t) = d T y = d^y = at r = 0, 
thus from Eq. (jl^) one can see that the sweep can be stopped at t = and the system will 
remain in this state. To demonstrate this behavior one should rewrite the wave function of 
the system tp — ci(t)-0i + c^it)^^ determined above in the basis of the exact states 



^±(t) = [±K±{t)ih + K T (t)^ 2 ] , K±{t) 



l± , W ® (26) 



\| sjl + w\t) 



corresponding to the time-dependent "energy eigenvalues" of Eq. @ as ijj = c + (t)ip+(t) 
c-(t)ip-(t). The corresponding probabilities are given by 



One can see that in resonance at t = (W — x = 0, cosy = ±1) the situation depends on the 
sign of cosy. For sweep functions given by Eq. jl^ ) with a < a max one has cosy = — 1 and the 
system is in the ground state (|c_(i)| 2 = 1, |c + (i)| 2 = 0). For the special sweep functions with 
a = a m ax one has cosy = 1 and the system is in the excited state (|c_(i)| 2 = 0, |c + (t)| 2 = 1), 
the transition from the ground state ip- (t) to the excited state ip+ (t) occuring for n ^> 1 in the 
vicinity of u — w ma x of Eq. ([la). These features are illustrated for the sweep functions of Eqs. 
( p^| ) and ( pl| ) in Figs, pjc andpji. Note that for n ~ 1 and a ~ a max = 1 the sweep function of 
Eq. ( pT| ) returns to — oo at r = oo and the system ends up in the excited state = ip2- 

To conclude, we have solved the inverse Landau-Zener problem for a two-level system, i.e., 
we have established the form of the sweep function W(t) = ei(t) — £2{t) needed to insure a 
given probability |ci(£)| 2 to be in the state 1 at time t. We have discussed two examples for 
which the system performs a complete transition from 1 to 2. In contrast to the direct LZ 
problem with a linear sweep, this can be accomplished by a finite sweeping rate. Further we 
have found such sweep functions that prepare the system in its exact ground state or in its 
exact excited state at resonance. It is not difficult to find W(t) that prepare the system in an 
arbitrary mixture of the exact states at resonance, which can be potentially used in quantum 
computing based on tunneling units. 
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